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Abstract. Tilted Bianchi type I cosmological models filled with disordered radia- 
tion in presence of a bulk viscous fluid and heat flow are investigated. The coefficient 
of bulk viscosity is assumed to be a power function of mass density. Some physical 
and geometric properties of the models are also discussed. 
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1. Introduction 

General Relativity describes the state in which radiation concentrates 
around a star. Klein (1948) worked on it and obtained an approxi- 
mate solution to Einsteinian field equations in spherical symmetry for 
a distribution of diffused radiation. Many other researchers (Singh and 
Abdussattar, 1973; Roy and Bali, 1977) have worked on this topic and 
obtained exact static spherically and cylindrically symmetric solutions 
of Einstein's field equations with exception as well. Roy and Singh 
(1977) have obtained a non-static plane symmetric spacetime filled with 
disordered radiation. Teixeira, Wolk and Som (1977) investigated a 
model filled with source free disordered distribution of electromagnetic 
radiation in Einstein's general relativity. 



The general dynamics of tilted models have been studied by King 
and Ellis (1973) and Ellis and King (1974). The cosmological models 
with heat flow have been also studied by Coley and Tupper (1983. 
1984); Roy and Banerjee (1988). Ellis and Baldwin (1984) have shown 
that we are likely to be living in a tilted universe and they have in- 
dicated how we may detect it. Beesham (1986) derived tilted Bianchi 
type V cosmological models in the scale-covariant theory. A tilted cold 
dark matter cosmological scenario has been discussed by Cen, Nickolay, 
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Kofman and Ostriker (1992). 

The majority of the studies in cosmology involve a perfect fluid. 
However, observed physical phenomena such as the large entropy per 
baryon and the remarkable degree of isotropy of the cosmic microwave 
background radiation suggests analysis of dissipative effects in cos- 
mology. Furthermore, there are several processes which are expected 
to give rise to viscous effects. These are the decoupling of neutrinos 
during the radiation era and the decoupling of radiation and matter 
during the recombination era. Bulk viscosity is associated with the 
GUT phase transition and string creation. The model studied by Mur- 
phy (1973) possessed an interesting feature in that the big bang type 
of singularity of infinite spacetimc curvature does not occur to be a 
finite past. However, the relationship assumed by Murphy between 
the viscosity coefficient and the matter density is not acceptable at 
large density. The effect of bulk viscosity on the cosmological evolution 
has been investigated by a number of authors in the framework of 
general theory of relativity (Padmanabhan and Chitre, 1987; Johri and 
Sudarshan, 1988; Maartens, 1995; Zimdahl, 1996; Pradhan, Sarayakar 
and Beesham, 1997; Kalyani and Singh (1997; Singh, Beesham and 
Mbokazi, 1998; Pradhan et al., 2001, 2002). This motivates to study 
cosmological bulk viscous fluid model. 

Recently Bali and Meena (2002) have investigated two tilted cos- 
mological models filled with disordered radiation of perfect fluid and 
heat flow. Meena and Bali (2002) have obtained two conformally fiat 
tilted Bianchi type V cosmological models. Very recently tilted Bianchi 
type I cosmological model for perfect fluid distribution in presence of 
magnetic field is investigated by Bali and Sharma (2003). In this paper, 
we propose to find tilted Bianchi type I cosmological models filled with 
disordered radiation in presence of a bulk viscous fluid and heat flow. 



2. Field Equations 

We consider the Bianchi type I metric in the form 

ds^ = -dt'^ + A^dx'^ + B'^dy'^ + C'^dz'^, (1) 

where A, B and C are function of t only. 
The Einstein's field equations are given by 

- = -SttT/ (c = 1, G = 1 in gravitational unit) (2) 
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where is the Ricci tensor; R = g^^Rij is the Ricci scalar; and T/ is 
the stress energy-tensor in the presence of bulk stress given by 

T/ = (p + p)viV^ + pgi + QiV^ + Viq> , (3) 

and 

P = P- (^) 
Here p, p, p and ^ are the energy density, isotropic pressure, effective 
pressure, flow vector, bulk viscous coefficient respectively and Vi is the 
flow vector satisfying the relations 

gijvV = -1 (5) 

QiQ^ > 0, (6) 

qy = 0, (7) 

where Qi is the heat conduction vector orthogonal to Vi. The fluid flow 
vector has the component ( ^^'^^ , 0,0, cosh A) satisfying Eq. (5) and A 
is the tilt angle. 

The Einstein's field equations (2) for the line element (1) has been set 
up as 

„ r/ N . , 9 » sinhA- Baa Caa BaCa , , 

-SttUp + p) sinh^ A + p + 2gi — ^] = -ii + + -i^, (8) 

_8;r^.= — + — + -^, (9) 

Aaa Baa AaBa , , 

-8.p=^ + ^ + -l^. (10) 

(p + p\A sinh A cosh A + gi cosh A + ^^^^ = 0, (12) 

cosh A 

where the suffix 4 at the symbols A^ B, C denotes ordinary differenti- 
ation with respect ti t. 



3. Solutions of the field equations 

Equations (8) - (12) with Eq. (4) are five independent equations in eight 
unknowns A, 5, C, p, p, ^, q and A. For the complete determinacy of 

the system, wc need three extra conditions. 

We assume the model is filled with disordered radiation which leads to 

p = 3p (13) 
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and 

A = {B C)" (14) 
where n is any positive real number. Eqs. (8) and (11) lead to 

-B44 , C44, 2B,iCi Aid . A.iBi 

Using Eq. (13) in Eq. (15) reduces to 

Baa Caa IBaCa AaCa AaBa , „. , . 

Equations (9) and (10) lead to 

.644 Aj Cj N _ p. 

c~b a^c~b'~ ^ ^ 



which leads to 



where BC = e, ^ = u and k is a constant of integration. 
Eqs. (9) and (10) also give 



(18) 



2^44 S44 C44 yl4C4 ^4^4 

^ + ^ + ^ + ^ + ^ = -167r(p-2C0) (19) 

From Eqs. (16) and (19), we obtain 

Aaa Baa Caa BaCa AaCa AaBa „ /„„^ 



which can be rewritten as 

(„ + i)!li + „4_i!| + i4 = o (21) 

e 4e2 4z/2 

where A = e". From Eqs. (18) and (21), we get 



ff^ I (22) 

4(n + l)e 4(n + l)e2"+i ^ ^ 

where 64 = /(e). Eq. (22) leads to 

n 4n+l 

2 _ + (4n + l)/cie^^ 
" (4n + l)e2n ^^^^ 

where ki ia a constant of integration. On integrating Eq. (18) we obtain 

log V = k^{An + l) [ == (24) 

eyfe^ + ^i(4n + l)e2(n+i) 
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Hence the metric (1) reduces to the form 



where v is determined by Eq. (24). 
By using the following transformation 

e = T, X = X, y = Y, z = Z 

the metric (25) takes the form 

(4n +1)T^'' 



(25) 



ds 



where 



F + (4n + l)kiT'^^^ 



1/ 



/o5( V = k^J{An + 1) / ^ 



4n+l 

A;2 + A;i(4n + l)r2("+i) 



The pressure and density of the model (26) are given by 

(4n + l)A;i 



(26) 

(27) 



87r(p-C^) = 

87r(p - 3^^) 

The tilt angle A is given by 

cosh^ A 



4n2+4n + 3 

8(n + l)T 2(n+i) 
3(4n + l)fci 

4n2+4n+3 

8(n + l)r 



2n + 3 



Hence 



sinh^ A 



2(2n+ 1) 

1 - 2n 
2(2n+ 1) 



(28) 

(29) 

(30) 
(31) 



The scalar of the expansion 9, calculated for the flow vector v*, is given 

by 



(n + l) 

J'n+l 



4n+l 



(2n + 3)[/c2 + (4n + l)A;iT5(^] 
2(2n+l)(4n + l) 



(32) 



If we put ^ = in Eqs. (28)- (29), we get the solutions as obtained by 
Bali and Meena (2002). 

Thus, given ^(t) we can solve the system for the physical quantities. 
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Therefore to apply the third condition, let us assume the following 
adhoc law (Maartens, 1995; Zimdahl, 1996) 



m = ^o/ 



(33) 



where and m are real constants. If m = 1, Eq. (33) may correspond 
to a radiative fluid (Weinberg, 1972), whereas m = | may correspond to 
a string-dominated universe. However, more realistic models (Santos, 
1985) are based on lying the regime < m < 5. 



3.1. Model I: 



(e = eo 



When m = 0, Eq. (33) reduces to ^ = ^0 and hence Eqs. (28) and (29) 
with the use of Eq. (32) load to 



_ (n + 1)^0 



(2n + 3)[fc^ + (4n + l)fciT^ 
2(2n + l)(4n+l) 



4n + l 



n+1) 



- + 



(4n + l)ki 



P 3rn+l '\ 



4n+l 



(2n + 3)[fc2 + (4n + l)kiT^(^] 



4n2+4re+3 

64(n + l)7rT 

(34) 

(4n + l)ki 



2(2n+ l)(4n + 1) 



4n2+4n+3 

192(n + l)7rr 2{"+i) 
(35) 



3.2. Model II: (^ = ^op) 

When TO = 0, Eq. (33) reduces to ^ = ^op and hence Eqs. (28) and (29) 
with the use of Eq. (32) lead to 



P 



(4n + l)A;i 



271^ + 1 



647r(n + l)T2("+i) 



T('^+i) -3(n + l)CoW(2n + 3) 



(4n + l)ki 



4n+l 

fc2 + (4n + l)/ciT2{"+i) 
(36) 



2n2 + l 

1927r(ra + l)r2(r.+i) 



r("+i)-3(n + l)CoW(2n + 3) 



4n+l 

/s2 + (4n + l)fcir2(n+i) 
(37) 
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4. Some Physical and Geometric Properties of the Models 

The flow vector and heat conduction vector for the models (26) 
are obtained by Bah and Meena (2002) 



V 



1 / (1 - 2n) 



Y 2(2n + 1) 



I 2n + 3 
2(2n + 1) 



6(4n + l)(2n + 3) / (1 - 2n) 



647r(n + 1)T^ 



-^^V2(2n+1)' 



_ 6(4n + 1)^(1 -2n) / 2n + 3 
~ 647r(n + 1) V 2(2n + 1) 



2n^+2n-3 
J" 2{n + l) 



4n + l 



a2 + 6(4n + 1)T^^ 



(38) 

(39) 
(40) 

(41) 



The non- vanishing components of shear tensor (ajj) and rotation tensor 
(iOij) are obtained as 



0-11 



0'22 = 1^1 



(2n - 1) 



■2n + 3 



2n + 1 



nn-l 



4n+l 

a2 + 6(4n + l)r2("+i) 



\ 2(4n + 1) 



(42) 



2n + 3 



, / 4n + l 

3aV4n + 1 + (1 - 2n)y + 6(4n + 1)T^(^ 

(43) 



2(2n + l)(4n + l 
where v is already given by (27). 

1 r 



2n + 3 



/ 4n + l , 

1 - 2n) + &(4n + 1)T^(^ - 3oV4n + 1 

(44) 



(744 



Cl4 



i^Y2(2n + l)(4n + i; 

(l-2n)2 / 2n + 3 



" 6(2n + 1) Y 2(2n + 1) (4ra + 1) 



+ 6(4n + l)r2{"+i) 



(6 - 2n - 32n^) / (1 - 2n) 
12(2n + l)i Y2(4n + 1) 



4n + l 



a;i4 = ni 



(1 - 2n) 



2(2n + l)(4n + 1) 



+ 6(4n + l)r2^^ 



a2 + 6(4n + 1)T^(^) 



(45) 
, (46) 

(47) 
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The rates of expansion Hi in the direction of X, Y, Z-axes are given 

by 

■ / 4TI+1 

+ 6(4n + 1)T^(^ 



n 



V4n + 1 



J'n+l 



Ho 



Ho 



+ 6(4n + 1)T^(^ + ay/ An + 1 
2V4n + ir"+i 



4?x+l 



a2 + 6(4n + 1)T^^^ - av'4n + 1 
2V4n + IT'^+i 



(48) 
(49) 

(50) 



The models in general represent shearing and rotating universes. The 

expansion in the models decreases as time increases and the expansion 
in the models stops at T = oo. There is a big bang in the models at 
r = Oifn + l>0. There is no rotation in the models for n = but 
it is shearing and goes on decreasing as time increases. For n = 0, the 
Hubble constant Hi = 0, H2 7^ 0, H^ 7^ 0. The fluid velocity vectors 



V 



^ and u*^ for n = 0, are given by f -"^ = 



V 



|. Both density and 
pressure in the models become zero at T = 00. Since limf^oo ^ 7^ 0, 
the models do not approach isotropy for large values of T for general 
value of n and n = also. 



5. Particular Models 

We consider the metric 

ds^ = -dt'^ + dx^ + B^dy'^ + C^dz'^, (51) 
where B, C are functions of t alone. 

Here we only assume that the model is filled with disordered radiation 
which leads Eq. (13). Following Bali and Meena (2002), we obtain, from 
Eqs. (51), (2) and (3) 

B^^tlif^, ,52) 

K2 



where k, k\ and ^2 are constants of integration. Here BC = e, ^ = u, 
€4 = /(e) and 

^ = (54) 



abhasty.tex; 7/02/2008; 2:16; p. 8 



TILTED BIANCHI TYPE I COSMOLOGICAL MODELS 



9 

(55) 



Thus, after the suitable transformation of coordinates, the metric (51) 
reduces to the form 

1 

ds'^ = (T^ - k'^fdT'^ + dX'^ + (T - kfdY'^ + (T + kfdZ'^ (56) 

Kg 



The pressure and density for the model (56) are given by 



^2 



647r(r2 - /fc2)2 



{p-m 

The tilt angle A is given by 



3^2 



647r(r2 - fe2)S 



(57) 

(58) 



coshA = y| (59) 
The scalar expansion calculated for the flow vector u*, is given by 

(60) 



fool 


klT 


V 2 


(r2 - fe2)3 



Thus, given ^{t) one can solve the system for the physical quantities. 
5.1. Model I: (^ = ^0) 

When m = 0, Eq. (33) reduces to ^ = ^0 and hence Eqs. (57) and (58) 
woth the use of Eq. (60) lead to 



P 



(r2 - A;2)2 



^2 



2 ^ 647r(r2 - A;2) 



>^ (T2 - A:2)2 
5.2. Model I: (^ = ^op) 



2 ^ 647r(r2 - A;2) 



(61) 

(62) 



When m = 1, Eq. (33) reduces to ^ = ^oP and hence Eqs. (57) and (58) 
with the help of Eq. (60) lead to 



ko 



P - 



327r(T2 - P) 2(r2 - A:2)2 - 3^6 k^T 



(63) 
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3A^9 



327r(T2 - ik2) 2(r2 - fe2)2 _ 3^ 



(64) 



6. Some Physical and Geometric Properties of Particular 

Models 

The non- vanishing components of shear tensor cr—ij and rotation tensor 
ujij for the models (56) are given by (Bah and Meena, 2002) 



1 /3 



^^^ = -2V2 



ki /3 



1 /3 



(T44 



<Ti4 



1 

2^2 



(T2 - P)2 

{T + kf 

{T-3k){T + kf 
(T - A;)2 

(T2 - A;2)2 

(r2 - A;2)2 



(65) 

(66) 

(67) 

(68) 

(69) 
(70) 



UJl4 = 

The expressions for vi, V4, qi and for the metric (56) are obtained 
as 

(71) 



Vl 



16 



V2' 



v,=-f,iT^-kr 

"-2 



3 
2' 



91 



94 



3^2 



64V2 7r(r2-a2)3^ 

3^2 



(72) 
(73) 
(74) 



64^6 ^(r2 - a2)3 • 
The rate of expansion Hi in the directions of X, Y, Z-axes is given by 

Hi = 0, (75) 



Ho 



^2 



2(r + A;)(T- A;)2' 



(76) 
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" 2(r- A;)(r + A;)2- ^'^^^ 

The models start with a big bang at T = ztk and the expansion in the 
models decreases as time increases and the expansion in the models 
stops at T = 0. The x-component of the Hubble parameter is zero due 
to the assumption of metric. However, the y and x components become 
infinite at T = A;. The models, in general, represent shearing and non- 
rotating universe. 
Since 

lim 5 7^ 0. 

t->oo 9 

Hence the models do not approach isotropy for large values of T. When 

T ^ k, qi ^ —oo and is constant. There is a singularity in the 
models at T = k. This singularity is pan cake type (MacCallum, 1971) 
as gii and 522 is zero at the singularity T = k. 
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